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ABSTRACT :.: 


Measurements of three wave surface profiles in a line 
perpendicular to the beach extending across the breaker line 
and a set of horizontal and vertical water particle velocities 
at the breaker line were taken in a mild surf and spilling 
type breakers. The application of a linear wave theory 
transfer function to the surface profile spectrum resulted 
in a calculation of the horizontal water particle velocity 
spectrum which was 13 percent lower than the measured 
spectrum peak and showed a phase five degrees ahead of the 


measured horizontal water particle velocity. Celerity 


frequency range of the spectra showed that the wave phase 
speed was not frequency dependent and was predictable by the 
linear theory shallow water approximation. The finding that 
a dispersive linear wave system predicts the wave kinematics, 
while a non-dispersive system was observed, suggested a 
contradiction and was discussed. Capacitance wave gages 

and electromagnetic flow meters were found to be excellent 


AnNStrUmMentsS for usemm the surf zene. 
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I. INTRODUCTION 


A. OBJECTIVES OF THE EXPERIMENT 

The kinematics of wave forms and water particle vel- 
ocities in the surf and breaker zones has been an enigma 
to science for many years. Laboratory studies have estab- 
lished the fact that the linear solution to the wave 
equations, which applies well to waves of moderate height 
in deep water, cannot fully describe monochromatic waves 
and associated water particle velocities during the break- 
ing of waves. This experiment was designed to measure 
surface profiles and water velocities in the breaker zone 
and compare the results with present theory. AS a secondary 
objective the accuracy and usefulness of capacitance wave 
gages and electromagnetic flow meters in the surf environ- 
ment were to be evaluated. The experiment also added to 
a data bank of measurements taken in the surf zone for 
studies being conducted at the Naval Postgraduate School. 

It is expected that the data obtained will ultimately 
contribute to the establishment of a theory which will 
accurately predict various wave parameters from one set of 


measurements taken in the breaker zone. 


B. APPROACH TO THE PROBLEM 
Measurements of water surface profiles at four locations 
in a line perpendicular to the beach and two sets of 


horizontal and vertical water particle velocities at the 
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breaker point were attempted during this research. An 
energy-density spectrum was catculadted for each data set 
and coherence and phase were calculated between data sets. 
A wave profile spectrum was converted to horizonal and 
Vertical velocity spectra using linear theory transfer 
functions. These spectra compared with the measured ve- 
locity spectra taken at the same point. Solitary wave 
theory was also compared to the measured EESULES,. \the 
assumption that the processes in the breaker zone are sta- 
tionary and ergotic was investigated. 


¢ 


C. PREVIOUS WORK 

This experiment is a continuation of work conducted at 
the Naval Postaraduate School by Steer (1972) ana 
Richardson (1973) under the direction of Thornton (1968, 
1972). The problem of taking direct measurements in the 
Surf zone was first attacked by Inman (1956) who infered 
water particle velocities from measured drag forces under 
the waves. Miller and Zeigler (1964) measured in situ 
particle motions using acoustic and eleeceromagnetle Current 
meters. Walker (1969) made the same type of measurements 
using the propeller type flow meters. Wood (1973) measured 
waves and currents in the surf zone using movies of dye 
movement and capacitance wave gages. Work by Dean (1974) 
gives an analysis of various wave theories applicable in the 


Surt zone. 
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II. DATA COLLECTION AND PREPARATION 
RR EARALLION 


A. EQUIPMENT 

Measurements were taken using two electromagnetic flow 
meters, three capacitance wave height gages, and one 
resistance wave height gage. All equipment was calibrated, 
in the laboratory, prior to the experiment. 

The two flow meters were Marsh-McBirney Model 72] 
Electromagnetic Current Meters. The flow meter operation 
is based on Faraday's principal of electromagnetic induc- 
tion. Each probe measures water velocity in two orthagonal 
directions, horizontally (shoreward) and Vertically: 
Measurement accuracy was determined to be +0.02 m/s during 
calibration. A carpenters level was used to establish axis 
alignment in the field with an estimated error of Bs 
degrees. 

The flow meters were calibrated in the laboratory using 
an oscillating platform attached to an GCCeEntric armus0.7 20cm 
long driven by a variable Speed motor using the method of 
Steer (1972) and Richardson (1973). The speed of the motor 
was determined using a stopwatch to measure the time re- 
quired for five revolutions and then making a comparison 
with the average of ten maxima voltage outputs measured on a 
Strip chart recorder. Because the platform became unstable 
at velocities above 1.3 m/s (0.4 Hz), calibration stopped 


at this point. The flow meter calibration plots are shown 
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in Figure (1). A linear regression was performed on the 
actual velocity versus voltage measurements to obtain the 
calibration factors listed in Table (1). 

The three capacitance wave staffs were constructed by 
the author using the design of McGoldrick (1969). The 
2-1/2 m gages consisted of a 1.0 cm diameter stripped 
coaxial cable mounted on a painted wooden frame as shown in 
Figure (2). The capacitance wave gage operates on the 
principal that a change of capacitor plate dimensions 
changes the circuitry voltage output. The insulated coaxial 
cable and the sea water act as the plates and the cable 
insulation as the dielectric. As the wave elevation changes, 
the Capacitance of the circuit changes and the voltage output 
responds linearly. 

The frame was rather heavily constructed and attached 
close to the tower to maintain a fairly rigid probe wire and 
to reduce water resistance on unattached parts. When assembled, 
the frame became somewhat unwieldy in the surf and created 
some resistance to flow around the support tower. The accuracy 
of the gages, however, was not believed effected. The 
electronic gear, including power source, was mounted ina 
waterproof brass container on top of the frame and two conductor 
cables carried the output signal shoreward. 

The wave gages were calibrated in the laboratory to a 
depth ot i 7 mw Accuracy was eStimated to be +0.005 m. The 


calibration plots are shown in Figure (1) and the results 
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TABLE (1) - INSTRUMENT DESIGNATIONS 
AND CALIBRATION INFORMATION 


= eT i i is Se re ey * . 


a 


TRACK # 





Outer 
Profile 


TLE 





Horiz/Vert 
Velocity 






INSTRUMENT FM1X/FM1Z 















MANUFACTURER Baylor Marsh- 
McBirney 
MODEL/SERIAL# = 13528R 217 95 










RANGES 
3.40-8.29v} 2.40-5. 35v| 0.33-0.16 +1.0v 
0.0-1.52m° +3.05m/sec!| 


VOLTAGE OUT 
MEASUREMENT 


HEIGHT ABOVE 
BOTTOM : , 0.24m Or om 











LAB CALCULATIONS! 
MULT FACTOR 
ADD FACTOR 

STD. ERROR 







-9.405 +10.908 
+3.078 -0.223 
+0.014m/sec!} 












-1.753 
+0.004m 






PROCESSING GAINS“ 





FINAL CALIBRATION 


+0.024 
-0.068 


MULT FACTOR i) 0G POL O7 2 
ADD FACTOR eS me 69 





1 Prom linear regression of laboratory calibrations as 
Shown in Figure (1). 


2 A - Instrument to analog tape, B - Analog tape, through 
filtering, to digital tape, and C - 1/(AxB), applied to 
digitally recorded voltage to compute instrument output. 


3 Pactors Apemedeto data in Gigital form (including 


instrument height) during computer analysis to obtain actual 
surface height and velocity values. 
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CAPACITANCE WAVE GAGE CONSTRUCTION 


ELECTRONICS IN SEALED 
BRASS CONTAINER 


DETAIL 1 


a) 
J HOLLOW BRASS 
INPUT A otup with 
WIRE SOLDERED 
DETAIL 2. DETAIL 1 INISIDE 
Cue d | 7 ’ 


NUT ANO 
Are WASHER TO 
TENSION 
CABLE 









CABLE 
CABLE 
(| 
eS 
DETAIL 2 
250 cm, 
WATE R PROOF 
| SEALER 
| ; 
1] LOOP OVER 
5 WOODEN 
Peg AT 
CABLE ENO 
lines eee 
a ee DETAIL % 
ol em, Jo 
TOWER FRAME - REDWOOD (2%4") SEALED 


WITH WATERPROOFING 


CABLE — CENTER OF RG -11 COAXIAL 
CABLE (APPROX. 1 CM OD) 


Paugure. (2) 
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of a depth versus voltage output linear regression are 
shown in the calibration factors of Table (1). 

The resistance wave gage was a Baylor Wave Measuring 
System Model #13528R with a 1.52 m wave height capacity. 
The gage consisted of two 1.27 cm stainless steel wire 
ropes mounted under tension on a tower and an electronics 
box located ashore. As the water level rose on the cables, 
they were shorted and a voltage was produced which was 
linearly proportional to the amount of exposed cable. 

Calibration consisted of simply shorting the cables at 
various intervals and recording the output voltages. Cali- 
bration accuracy was +0.02 m when measurements were taken 
away from either end of the gage. The calibration plot is 
shown in Figure (1) and the calibration factors determined 
by linear regression are shown in Table (1). 

All data was recorded on a 14 channel Sangamo Analog 
Recorder. Voltage inputs were reduced to a range of -1.0 
to +1.0v and recorded continuously. A strip chart recorder 
ys used to monitor both input and recorder output voltages 
during the experiment. One wave profile data set (WS3) was 


lost during recording. 


B. BEACH SET UP 

The experiment was conducted at Del Monte Beach at the 
Naval Postgraduate School Beach Laboratory on Thursday, 
23 May 1974. Data was recorded from 12:50 to 14:30 PM (PDT) 


(100 minutes). High tide was at 13:40 PM. The weather was 
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clear, with some gusty onshore wind (less than 10 knots) and 
no visible sea. The swell surf was mild with observed 
significant wave heights of between 0.5 and 0.75 m and periods 
of approximately 8 to 10 seconds. Breakers were predominantly 
spilling with an occasional plunging type. Refraction along 
the shore consistantly caused waves to break paralled to the 
beach. 

The four towers used were steel pipes 6.3 cm in diameter 
and 3.6 m high with a 1.0 m diameter base. They were held up 
with 0.9 cm polypropylene rope tied to homemade screw-type 
anchors. A typical tower arrangement is shown in Figure (3). 
The towers were placed on a line perpendicular to the shore 
at specified distances. The towers remained visibly stable 
during measurements and it was felt that tower movement was 
not appreciable. The polypropylene lines had a tendancy to 
stretch and loosen and knots did not hold well. The use of 
steel guy wires would have reduced line slackness problems. 
The anchors were difficult to both install and remove and 
did not hold unless deeper than about 0.5 m. The experiment 
was concluded when Tower Two became unstable and fell as a 
knot failed. Flow meter number two (FM2) was mounted too 
high on the tower and was occasionally exposed, causing 
false readings. Thus this data could not be used for ae 
continuous analysis. Figure (4) is a diagram of the beach 


profile and the instrument arrangement. 
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C. DATA PREPARATION 

Data was recorded at the beach on eight channels ata 
speed of 1 7/8 in/sec. The data was played back at 
7 1/2 in/sec and sampled at 100 Hz (1 sample per 0.01 
seconds) or 25 Hz real time (1 sample per 0.04 seconds). 
The analog to digital conversion was performed on a CDC/9300 
computer. A summary of tape recorder and digitization gains 
is included in Table (1). Playback data was amplified 
once, filtered, amplified again, digitized, and recorded. 
A matched set of Krohn-Heit filters with a cutoff frequency 
of 40 Hz (10 Hz real time) was used to low pass filter the 
analog signals. Filtering was conducted during this phase 
to reduce the possibility of 60Hz aliasing. The real time 
sampling interval was 0.04 seconds so the recorded Niguist 
frequency was 12.5 Hz. Digitized data was recorded on a 
nine track magnetic tape (labeled BUB191) in one file of 
597 records, each record consisting of 2048 points or 256 
data points for each of eight channels. Figure (5) gives 


examples of the recorded data taken during the experiment. 


PAA 





pepryH (m) 


DEPTH (m) DEPTH (A) 


VELOCITY (M/S) 


VELOCITY (™/8) 


SAMPLES OF DATA AS RECORDED 


4 
3 
; 4 OUTER PROFILE 
(WS2) 
4 
r) 
Ke) 
2 2.CENTER PROFILE 
4 (WS1) 
‘e) 
, 3. INNER PROFILE 
(BAY) 
(ee Nerf DONA... 
O 
3 
2 
$ 4, HORIZONTAL 
o NECOCTVY 
er (FMLYX) 
-2 
->3 
3 
: & VERTICAL. 
0 VELOCITY 
Et (FM1t) 
-2 
3 


Oo 5S» Jo 20 30 AO te) 
TIME SCALE (SECONDS) 


Figure (5) 
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Iifl. SPECTRAL ANALYSIS 


A. APPROACH TO DATA ANALYSIS 

Continuous time series records of all measurements were 
taken for 100 minutes. A thirty minute sample record of 
each data set was spectrally analyzed on an IBM/360 
computer. The power spectra generated described the com- 
position of the data over a frequency range in terms of its 
mean square value. An examination of a power spectrum 
indicates the regions of greatest mean wave or velocity 
amplitude and their corresponding frequencies and will give 
a general relationship of energy concentrations across the 
frequency range. A comparison of two spectra by means of 
a cross-spectrum and subsequently a coherence and phase 
plot with frequency will indicate the regions and degree of 
linear relationship and phase between two data sets. Co- 
herence can be thoughtof as the ratio between two measure- 
ments of the square of the gain factors for a constant 
parameter linear system. A coherence near 1.0 indicates 
that both spectra operate under the same gain factor or 
transform function, that is, their relationship is highly 
linear. A coherence near 0.0 indicates the two examined 
spectra are not related or are related only by a non-linear 
gain factor. Phase angle measurements indicate the time 


delay between two data set spectra. 
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The use of spectral analysis in the examination of data 
required the following assumptions be made: 

1. The data was random. In other words, the data used 
represented one of many possible results that could have 
occured or was one physical realization of the random data 
process. | 

2. The random process was stationary. That is, the 
data set could be described in terms of certain average 
values which are true for all data sets of the same process. 

3. The stationary random process was ergotic or all 
properties of the process for all time could be described 
by performing analysis on one finite data set. 

Many processes found in nature can be treated as being 
ergotic, stationary, and random. The validity of this 
assumption as it applied to this experiment will be examined 


in Section (1IV,C) of this paper. 


B. COMPUTER ANALYSIS METHOD 

A complete flow chart of data acquisition and analysis 
is presented in Figure (6) and is discussed briefly below. 
The equations uSed are given in Bendat and Piersol (1971). 
The parameters used for the spectral analysis process were 
chosen to balance the longest record which the computer 
could reasonably analyze with the best resolution over the 
frequency range of interest. The frequency band of 0.0 to 
1.5 Hz was chosen to cover the region of gravity waves 


(0.033 to 1.0 Hz). Low pass filtering was used twice to 
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DATA ACQUISITION AND ANALYSIS FLOW CHART 
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Figure (6) 
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Minimize 60 Hz aliaSing. Folding of 60 Hz noise around the 
Niguist frequency would have appeared as peaks in the 
spectra at 0.625 and 0.9375 Hz. No such peaks showed up 
during analysis, demonstrating the success of the filtration 
employed. 

The same 30-minutes of raw data was analyzed for each 
of the five channels. During each program, two channels 
were chosen, digitally filtered, and decimated by eight. The 
data was low pass filtered using a 25 weight filter with a 
cutoff frequency of 1.2 Hz and terminal frequency of 1.8 Hz. 
This filter, illustrated in Figure (7), is given in Davidson 
(1970). The final analyzed time step was thus 0.32 seconds 
with a Niquist freguency of 1.5625 Hz. 

Hie Gata pCints were them converted to wave height 
and/or velocity values uSing the calibration factors of 
Table (1). A mean value was calculated for both data sets 
and the data was linearly detrended. The variance, standard 
deviation, and average period were, calculated. The average 
period was determined by finding the time between zero axis 
upcrossings; the calculated average periods appear low 
using this method. The mean values for each data set are 
listed in Table (2). 

For pairs of measured parameters, the auto-covariance 
functions were calculated using a five percent lag. Each 
auto-covariance was smoothed with a Parzen window. A 
Fourier transform was then applied to the smoothed auto- 


covariance functions and the two power density spectra 
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TABLE (2) - STATISTICAL PARAMETERS 


CALCULATED DURING ANALYSIS 


mRACK 


TITLE 


DESIGNATION WS2 WS1l BAY FM1X FM1X 


MEAN DEPTH 
OR VELOCITY 1.982 1.468 Lea23 +0.018 +0.048 
(m or m/sec) 


VARIANCE 













(m2 or m/sec?) 









STANDARD 
DEVIATION 


(m or m/sec) 
AVERAGE PERIOD Ses ee: 5.29 4.20 6.20 3.10 


(seconds) 
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determined. Bandwidth resolution for the spectra was 
0.00544 Hz. 

A cross-covariance function between data sets was then 
computed, smoothed with a Parzen window, a Fourier trans- 
form applied, and the cross-spectrum computed. The 
coherence and phase relations were calculated and the two 


spectra, the coherence, and the phase plotted. 
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IV. RESULTS 


A. INTRODUCTION 

The spectra generated were examined, the observed low 
frequency peaks were related to phenomenon which one would 
expect to find in the surf zone. The stationarity of the 
data sets was examined in order to determine the scope of 
applicability of the data. 

An examination of theories which may be used in the 
surf zone was conducted using two approaches. First, 
measured velocity spectra were compared with velocity spectra 
calculated from the center surface profile uSing a linear 
wave theory transfer function. Solitary wave theory was 
then to be used to develop a discrete line spectrum for a 
Similar type of comparison. However, the results showed 
that: (1) the linear theory transfer function calculated hori- 
zontal velocity very well and (2) there were no identifiable 
harmonics in the measured spectra which could be related 
to such a discrete spectrum; this indicated that solitary 
wave theory was not applicable under the experiment 
conditions. 

A second approach to the examination of applicable 
wave theories was to compare measured values of celerity or 


wave phase speed with theoretical values. 
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B. EXAMINATION OF LOW FREQUENCY SPECTRAL COMPONENTS 

Figures (8) and (9) give the results of the spectral 
analysis when comparing the outer and center profiles and 
the center and inner profiles. Figure (10) is a spectral 
comparison between horizontal and vertical measured ve- 
locity spectra. Figure (11) is a superposition of the three 
surface profile spectra. Table (3) is a summary of low 
frequency spectral peaks found throughout the five spectra 
generated and gives a wave height or velocity value for 
each peak relative to a maximum spectral peak of 1.0. These 
values are presented in this manner for comparison purposes 
only. 

The following observations can be made concerning the 
results shown in the Figures and Table (3): 

1. The maximum wave energy for all spectra iS concen- 
trated at a primary frequency of 0.125 Hz corresponding 
to a period of 8.0 seconds. 

2. A secondary peak exists at 0.098 Hz (10.2 seconds). 

3. An approximate subharmonic of the primary frequency 
exists at 0.060 Hz (16.7 seconds). 

4. A fourth peak of interest shows up at 0.011 Hz 
(92 seconds) on the inner profile and vertical velocity 
Spectra. 

5. A major peak on the horizontal velocity spectrum 


only exists at 0.016 Hz (61.2 seconds). 
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SPECTRA - OUTER SURFACE PROFILE VS. CENTER SURFACE 
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SPECTRA ~- CENTER SURFACE PROFILE VS. INNER SURFACE. 
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SPECTRA - HORIZONTAL PARTICLE VELOCITY VERSUS 
VERTICAL PARTICLE VELOCITY 
AT 0.93 METERS 
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SUPERPOSITION OF THE THREE SURFACE PROFILE 
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6. From a qualitative examination of the high fre- 
quency ends of the three profile spectra (Figure (11)), it 
appears that turbulence increases to a peak at the breakers 
and the waves seem to reorganize slightly inside this point. 

Without further data analysis, discussion of these ob- 
servations can be made only in a general sense. The 
observation of a primary frequency at 0.125 Hz agrees well 
with the wave periods of eight to ten seconds observed 
during the beach experiment. The additional peak at 0.098 
Hz indicates that incoming waves were probably composed of a 
combination of 8.0 and 10.2 second waves. 

Guzza and Davis (1974) state that in a theoretical 
analysis of edge waves excited by incoming waves, the prom- 
inent edge wave mode is the first subharmonic of the primary 
frequency. An energy peak was observed at this frequency. 
The author feels that the peak at 0.011 Hz, which shows on 
the inner profile and vertical velocity spectra only, may 
also be caused by an edge wave in the swash zone. Further 
experiments with profile measurements in a line parallel to 
the beach will clarify edge wave observations. 

The comment that the turbulence inside the breaker zone 
decreases must be tempered by the fact that the inner wave 
gage was not the same type as the other two gages. It is 
‘felt that the Baylor wave gage was not as responsive to high 
frequencies as were the capacitance wave gages. Future 
studies with matched gages across the entire surf zone will 


clarify this question of wave organization. 
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C. EXAMINATION OF STATIONARITY 

As the tide rises and the mean water depth increases, 
the breaker zone should move shoreward, and the mean para- 
meters measured (mean depth and variance) at the towers 
should change. This series of events would result ina 
non-stationary situation. In Section (III.A) it was stated 
that if the data were ergotic, the measurements taken and 
analyzed during this experiment were typical of all waves 
in the breaker zone. This rather bold assumption can be 
examined by checking for stationarity of the data over the 
period of measurements. 

The data is assumed to be weakly stationary if the mean 
and mean square values calculated are constant with time or 
if the auto-correlation function or its Fourier transform, 
the energy density spectrum, is constant for adjacent 
lengths of records. Both approaches were taken to verify 
stationarity. Energy density spectra were computed for the 
outer and center profiles during two independent 30 minute 
time segments and are illustrated in Figure (12). Ina 
general sense (examining slopes and energy peaks), they 
compare favorably. 

Mean depths and variances were then calculated for ten 
overlapping samples of 30 minutes each and are plotted in 
Figure (13). A non-parametric run test was applied to the 
data as described by Bendat and Piersol (1971) to check for 
stationarity. On the basis of this test, the data did not 


qualify as stationary. 
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COMPARISON OF TWO !NOEPENDENT TIME SEGMENT 
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The non-stationarity can be explained by a number of 
Observations. First, it must be re-emphasized that the 
waves measured were unusually small and consistantly of 
the spilling type. The bottom hump shown in the beach 
profile (Figure (4)) at -18m appeared to trip the breakers 
in the vicinity of the center wave gage despite an apparent 
increase in the still water depth. The height of waves 
outside the breaker point, which can be related to the 
variance of the outer profile measurements, decreased as 
the depth increased with time. This was an indication of 
a decrease in the shoaling effect on the waves. Set down, 
the difference between still water level and mean depth at 
the breakers, increased as the breaker wave height decreased 
at the center wave gage. Future analysis should concentrate 
on this set down effect and its relation to set down theory. 

Thus, the rising of the tide was the main factor contri- 
buting to the non-stationarity of the data. Despite the 
fact that the data was non-stationary, the analysis was still 
valid when the fact that the data applies to a very special 


set of circumstances is considered. 
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D. LINEAR THEORY APPLICATIONS 


Airy or linear wave theory assumes that the surface 


profile is sinusoidal, l.e., 

nice Becos (ot) 
Linear theory horizontal and vertical 
velocities at the depth z (z positive 


mean water level are given by 


cosh (k(h+z)) 


ee) eee Sinh (kh) 


COS 


and 


_ Sinh (k(h+z)) -: 
w(t) =e Ao —Sinh (kh) sin 


water particle 


upward) 


(Ot) 


(ot) 


below the 


Horizontal and vertical velocity spectra can be computed 


from the wave spectrum using the equations 


c ea) = | S958 ee S 


uc sinh (kh) 

and 
s (f) = [sinh (kth+z))] ° 
we sinh (kh) 


The wave number, k, was calculated using an iteration of 


S (4) 


Srilee) 


the linear theory characteristic equation 


o¢ = gk tanh (kh) 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 





Calculated horizontal (Sia) and vertical (S,,.) velocity 
Spectra were computed using the center surface profile 
Spectrum and compared with measured spectra (S,,, and S,p) 
as shown in Figures (14) and (15). The calculated and 
measured horizontal velocity spectra, Syo and Sym, (Figure 
(14)) compare well, indicating that linear wave theory was 
applicable for this experiment. On the other hand, the 
calculated and measured vertical velocity spectra (S,,, and 
Suan) shown in Figure (15) did not agree well at all. The 
lack of agreement between these vertical velocity spectral 
magnitudes could be partially attributed to a low signal to 
noise ratio in the measurements (note the vertical velocity 
time sample relative to horizontal velocity time sample in 
Figure (5) which greatiy reduced the accuracy of the measured 
velocity spectrum. 

Phase measurements are used to substantiate the findings 
that linear theory was applicable. Equations (1) and (2) show 
that the center surface profile and. horizontal velocity 
Measurements should have been in phase (both have cos (ot) 
terms). An approximately zero degree phase was observed well 
beyond the 0.5 coherence point (Figure (14)). The cos (ot) 
and sin (ot) terms of equations (1) and (3) indicate that the 
surface profile and vertical velocity measurements should be 
90 degrees out of phase. Figure (15) illustrates that this 
phase difference was observed in a range up to the 0.5 


coherence. 
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SPECTRA — MEASURED VERTICAL VELOCITY VERSUS 
WERTICAL VELOCITY CALCULATED 
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Figure (15) 


45 





Finally, a high coherence would indicate a high degree 
of linearity in the transfer functions between different 
Spectra. Coherence above 0.5 was BESeHuad for an approxi- 
mate frequency range of 0.0 to 0.65 Hz (Figures (14) and 
(15)). This high value again substantiated the excellent 
results obtained when a linear wave theory transfer 
function was applied to the center surface profile spectrum 
when calculating the horizontal velocity. Whereas coherence 
is high for the same range, linear wave theory grossly 
overpredicts the vertical velocity spectrum. The author 
is led to believe that a linear relationship between the 
surface profile and vertical velocity exists but is not 
predicted well by linear wave theory. A further indication 
of the linearity between measurements was noted due to the 
high coherence between horizontal and vertical velocity 


measurements as shown in Figure (10). 


E. CELERITY COMPARISONS 

A summary of the various celerity calculations made in 
this section is included in Table (4). The measured 
celerity, Cis is computed considering the shift in phase 
existing between two adjacent wave gages as plotted in 


Figure (16). 
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At x) = 0, let 


ie Apeos eect ) (7) 
Then at Xo = xX 
No = Ag cos (kx-Ot) (8) 


Therefore, the phase shift, $6, between stations 1 and 2 


becomes 

a= Ee) x (9) 
By definition 

C(f) = Ey (10) 
Or substituting and solving for C, (£) 

C, (£) = ary & (11) 


The dotted lines in Figure (16) plot a constant change in 
phase with frequency over the coherent ranges of Figures 
(8) and (9). This surprising result indicates that wave 
celerity in the surf zone did not depend on frequency, that 
is, the waves were non-dispersive. 

The celerity was then computed for points between 
stations using the linear theory approximation for shallow 


water waves 


a 
Cy = (gh) * = constant (12) 
These values agree well with the measured celerity, C1 
(+ 4 percent). Again, as stated in the previous section, it 


appears linear theory is applicable in the surf zone. 
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Using the significant wave height, H,, calculated from 
the variance, celerity was computed using the approximation 
for both cnoidal theory and solitary wave theory 

C, = (g(htH,))* (13) 
These values overpredict the measured celerity by 20 percent. 

As a final step, celerity was calculated as a function 
of frequency using equation (10) for linear wave theory. 
The wave number, k, was calculated from an iteration of the 
linear theory characteristic equation, equation (6). The 
resulting celerities, C,(f), are graphed in Figure (17). 
The decrease in celerity with increasing frequency, indica- 
tive of a dispersive wave system and expected with linear 
theory, would result in a corresponding increase in phase 
shift with frequency. This relationship was not observed; 
rather, the phase shift with frequency was constant as 
illustrated in Figure (16). The value for C, (f) at the 
primary frequency, 0.125 Hz, agreed well with the measured 
celerity, C, (+ 2 percent). 

In summary, celerity calculations created a dilemma 
which future studies in the surf zone will have to resolve. 
The constant phase shift and measured celerity indicated 
that all frequency components are non-dispersive and ride 
on the back of a non-linear wave form moving at the speed 
C,- However, linear wave theory predicted the horizontal 
velocity and both C5 and Cqg (0.125Hz) closely approximate 
C,. If linear theory is applicable, the wave system is 


dispersive or each component moves at a celerity which 
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depends on its frequency. It appears that while linear 
theory can predict the relationships between waves it cannot 
be used to define the breaker zone wave form which is non- 


linear and non-disperSive. 


F. COMMENTS ON INSTRUMENTS AND BEACH EXPERIMENT 

As a corollary goal for this experiment, the instrumen- 
tation used and its placement in the surf zone was to be 
tested. The capacitance wave gages used were found to be 
reliable, inexpensive, hardy, and easily monitored in the 
surf zone. They were very accurate and appeared to give an 
excellent representation of wave profiles as they actually 
occurred. The work required to construct, calibrate, and 
install the gages was minimal. By mounting all electronics 
in a sealed container and running signal cables overhead, the 
chance of shorting the system and introducing noise was 
greatly reduced. 

The use of electromagnetic flow meters in the surf zone 
gave excellent results. The lack of moving parts which clog 
with sand and the ability to record and monitor signals on 
the beach virtually eliminated the problem encountered by 
previous experimenters of incomplete or questionable records. 
Calibration of the flow meters was somewhat difficult and 
time consuming and required special equipment. A new cali- 
bration scheme should be devised to increase the calibration 


frequency range and accuracy. The results obtained were 
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sufficient for the experiment, however. Calibration of 
the flow meters is required as their response, although 
linear, deviated from the manufacturers specifications, 
showing both a slight DC bias and a different calibration 
slope. 

The Baylor resistance wave gage also proved reliable, 
accurate, and very easy to calibrate. On the beach it was 
unwieldy and required a special tower to maintain tension on 
the cables. It was felt that the Baylor gage output was not 
aS accurate as that obtained with the capacitance gages. 

The selection of instrument locations on the beach was 
conducted four days prior to the actual meaSurements and 
was done by considering the beach profile and the location 
of breakers on that day. Wave gages were placed to span 
the breakers with the velocity measurements to be taken in 
the breaking waves. Fortunately, wave characteristics 
remained constant until the date of the experiment and 
placement was correct for the intended measurements. 

The key word for the experiment setup was organization. 
The equipment and manpower requirements for an examination 
of surf zone kinematics using temporary stations can 
quickly become overwhelming. It is felt by the author that 
five to six measurements are maximum for a one-man experiment 


iiweeieo Sure zone. 
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V. CONCLUSIONS 


With the conditions of a mild surf and spilling breakers, 
linear wave theory was shown to predict wave kinematics well. 
Spectral analysis showed a good agreement between the 
measured horizontal water particle velocity spectrum and 
the horizontal velocity spectrum calculated uSing a linear 
theory conversion factor on the surface profile spectrum 
over the highly coherent band of prominent wave energy. A 
similar high coherence and the agreement with predicted 
phase relationships over the same band found in a cross- 
spectral analysis of the surface profile and vertical water 
particle velocity led the author to believe that a linear 
relationship, although not predictable by linear wave 
theory, also existed which relates the vertical velocity 
spectrum to the surface profile spectrum. 

The shallow water approximation for celerity used in 
linear theory was an excellent prediction for wave phase 
speed in the experiment environment. A constant phase 
shift with frequency indicated the breaking waves were 
non-dispersSive or that astlestisy was not frequency dependent. 
Apparently all energy rode on the back of a non-dispersive, 
and therefore non-linear, waveform which was at the primary 
frequency. An inconsistancy thus arose in that linear theory, 
a dispersive theory, could be used to predict wave kine- 


matics; yet, celerity calculations all pointed to a 
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non-dispersive wave form. Future studies in the breaker 
zone are needed to resolve the question of linearity. 

The use of capacitance wave gages and electromagnetic 
flow meters is a convenient method of gathering accurate, 


continuous data in the surf zone. 
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